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COMPARISONS OF EXPERIMENTAL HYDROXYL RADICAL PROFILES 
WITH KINETIC CALCULATIONS IN A SUPERSONIC NOZZLE 


by Erwin A. Lezberg, Charles M. Rose, and Robert Friedman 
Lewis Research Center 

SUMMARY 

The hydroxyl (OH) mole fraction was determined from spectral line absorption mea- 
surements made in a supersonic nozzle that expanded the combustion products of hydro- 
gen and preheated air. Absorption measurements were made as a function of fuel- air 
equivalence ratio (0. 6 to 1. 2) at combustor pressures of 3. 6 and 1. 8 atmospheres at low 
supersonic area ratios (1.6 to 3. 55).) 

The measured OH mole fractions and previously reported static temperature mea- 
surements were compared to the output of a one-dimensional computer program that in- 
cluded finite rate chemistry. The conclusions drawn from these comparisons are as 
follows. 

The controlling reaction is 


H + OH + M * H 2 0 + M 

Increasing or decreasing the rate of this reaction by a factor of three from the "stan- 
dard rate” used in the calculations gives results compatible with the OH concentration 
measurements within the uncertainties involved in this measurement. Temperatures 
correlate best when a rate is used that is one-third the standard rate. 

The recombination reaction 


h + h + m*h 2 + m 

becomes important only when its rate is comparable to or larger than the rate of the OH 
recombination reaction. 

For lean mixtures a mechanism involving the reaction 

H + 0 2 + M * H0 2 +*M 

along with appropriate bimolecular HOg removal steps may be of importance since it of- 
fers an alternative recombination path. • 



INTRODUCTION 


The calculation of nonequilibrium nozzle expansions has received a great deal of at- 
tention recently, and a number of computer programs have been written for calculating 
performance. Reaction rate data needed for input to these programs have generally been 
obtained from shock tube measurements or experiments using flat flame burners. Mea- 
surements of atom or radical recombination rates for the hydrogen- oxygen reactions 
have generally been inferred in shock tube experiments from the reverse dissociation 
rates at conditions far removed from equilibrium. The application of these rates to a 
nozzle expansion remains somewhat in question for critical comparisons since spread in 
recent data of at least an order of magnitude exists for the important three-body recom- 
bination reactions (refs. 1 to 3): 


h + h + m*h 2 +m 

OH + H + M = H 2 0 + M 

(These are reactions (IV) and (V), respectively, as given in table I. ) 

The effects of third-body efficiency, particularly of water, are not well differenti- 
ated although some attempts at unravelling them in flat flame experiments have been 
made (ref. 3). 

Some nozzle experiments have been reported in which measurements of static pres- 
sures (refs. 1 and 4) or temperatures (refs. 5 to 7) are compared with exact or approxi- 
mate chemical kinetic calculations of performance. Reaction rates giving the best fit 
with the data have been inferred. In these comparisons, temperature measurements 
provide greater sensitivity to changes in rates than pressures. However, this does not 
allow separation of the effects of individual rates unless a single reaction such as (IV) 
or (V) is assumed to be controlling. 

With the addition of a second measured variable, such as a species concentration 
that is reasonably sensitive to the rates, it should be possible to test the sensitivity of 
changes in the temperature and in the species concentration to changes in the rate con- 
stants. Hopefully, some improvement in knowledge of the appropriate rates could then 
be achieved by comparisons with the measurements. 

The hydroxyl radical (OH) is known to be important to the kinetics of the hydrogen- 
oxygen reactions and is also spectroscopically accessible. Concentration of OH can be 
determined by spectral line absorption in the ultraviolet ^S + — electronic transition 
(ref. 8) by using known values for the transition probabilities. 
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The purpose of this investigation was to compare measurements of OH concentration 
and previously reported temperature measurements in a'nozzle (ref. 5) with machine 
calculations by using an exact kinetic calculation. Calculations were also used in estab- 
lishing the sensitivity of the calculated temperature and OH concentration to the various 
reactions. 

Combustion products of hydrogen and preheated air were expanded through a super- 
sonic nozzle. At a combustor pressure of 3. 6 atmospheres, inlet air temperature was 
in the range 1710° to 1915° K. At a combustor pressure of 1. 8 atmospheres, inlet air 
temperature was in the range 1725° to 1865° K. The range of equivalence ratios used 
was 0. 6 to 1. 2. 


EXPERIMENTAL TECHNIQUES 
Facility and Instrumentation 

A description of the experimental facility and basic instrumentation is given in ref- 
erence 9. Combustor inlet conditions for a Mach 6, subsonic-burning ramjet were simu- 
lated by using an alumina pebble bed to supply air at temperatures up to 1900° K and 
pressures up to 3.6 atmospheres to a water-cooled combustion chamber and nozzle. 

Hydrogen fuel was supplied from a cylinder trailer. It was injected into the combus- 
tion chamber through water-cooled fuel injection tubes drilled with 148 orifices 0. 040 
inch in diameter. The orifices were arranged to cover approximately equal areas of the 
combustor cross section. A DeLaval nozzle with a 3. 2- inch-diameter throat and a 7° 
half-angle divergent cone was equipped with 1/2-inch-diameter optical ports down- 
stream of the nozzle throat. The ports were purged with air during a run to keep the 
quartz windows clean. Static-pressure taps were drilled normal to the nozzle wall, and 
wall thermocouples were mounted flush with the inside wall. 


Hydroxyl Absorption Measurements 

The concentration of the OH radical was determined by line absorption in the ultra- 
2 + 2 

violet S — II electronic band (ref. 8). The description of the apparatus and method 
is given in some detail in reference 10. The apparatus is shown schematically in fig- 
ure 1. The light source is an end-view capillary discharge tube that is operated with a 
water vapor pressure of about 1. 0 millimeter of mercury. The emission line of the 
cooled source is absorbed by the broader absorption line of the combustion gases so that 
ideally only absorption at the line center is measured. 
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Monochromator 


Procedure and Data Reduction 


Transmission measurements were 
made in the 0-0 vibrational band by us- 
ing Q branch rotational lines from 
K = 12 to 19. This selection of rota- 
tional lines was made so that the absorp- 
tion followed the Beers’ Law Relation. 

The absorption measurements were 
made by first scanning past a spectral 
line of the lamp with a nonabsorbing gas 
(air) flowing through the nozzle and then 
by scanning the same line with the absorbing gas (hydrogen-air combustion products) in 
the light path. Since the light was chopped on the source side, and only the chopped sig- 
nal was amplified, possible emission from the gas was not recorded. 

The absorption at the line center is related to the ground electronic- state population 
by the f numbers for the transition and the shape of the absorption line. The f num- 
ber is defined as the ratio of the number of dispersion electrons to the number of ab- 
sorbers and is related to the electronic transition probability. (All symbols are defined 
in appendix A. ) 

The concentration of OH in the ground state is given by (ref. 10) 



5^d!1 
fa k t j'j" 


(hcto 


K 


k T 


x 2. 40x10 


12 


( 1 ) 


where y is the rotation-vibration partition function, b^ the Doppler half-width of 
the absorption line, and Aj^TjIjm the relative transition probabilities including the ef- 
fects of vibration- rotation interaction (ref. 11). 

The exponential term, which is the Boltzmann factor for the rotational distribution, 
was evaluated at the measured sodium spectral line (ref. 9) reversal temperature. The 
value of F, where 


F = f 


2J" 

A 


+ 1 
K 


-4 

is taken as 2. 0x10 ±10 to 15 percent (refs. 12 and 13). The absorption coefficient P 

is determined as the natural log of the measured reciprocal transmittance Iq/I, where 
Iq is the source intensity and I is the source intensity attenuated by the absorbing gas 
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The calculation of N QH from equation (1) is exact when assuming an infinitely narrow 
source width and pure Doppler broadening of the absorber line. Corrections for an as- 
sumed finite width of emission line and a combination of Doppler and collision broadening 
of absorber line can be made by replacing P‘ in equation (1) with the equivalent absorp- 
tion coefficient at the center of a Doppler broadened line where the integral of the absorp- 
tion coefficient over the Doppler line is equal to the integral for the actual line shape: 

p w 0" (2) 

D line 

The corrections are discussed in appendix B. The total correction increased the value of 
the absorption coefficient by about 20 percent at a chamber pressure of 3. 6 atmospheres 
and 17 percent at a pressure of 1. 8 atmospheres. 


Sources of Error 

Checks of rotational equilibrium were made by scanning a portion of the spectrum 
and determining the absorption coefficient 


for a number of lines in the 0-0 band. The rotational temperature can be determined 
by plotting log P^ y/ A K T j'j" against the rotational energy of the transition and by 

taking the slope of the resulting straight line. Within experimental error, the lines were 
straight and the temperature was generally within 100° K of the sodium reversal temper- 
ature measured at the nozzle centerline (ref. 9). The determination of the rotational 
temperature involved absorption measurements of some stronger lines that were 90 per- 
cent or more absorbed and, therefore, subject to errors because of self- absorption. 
Changes in Iq could occur through slight dirtying of the test section windows, through 
changes in source current or pressure, or through the movement of the source image 
with respect to the entrance slit due to vibration or refraction of the light by density gra- 
dients in the test gas. The change in the source intensity from these effects was usually 
less than 5 percent. The low-frequency noise level of the source was about 3 percent. 
Errors in the sodium reversal-temperature measurement could result in concentration 
errors due to the temperature effect on the population distribution, the half-width, and 
the gas density. A 5-percent error in the reversal temperature would result in a 4- 
percent error in the OH mole fraction for the Q^15 line at 1850° K. 
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Additional errors are introduced in the calculation of OH mole fraction by the uncer- 
tainty in the F value (10 to 15 percent), and in the corrections for the line shape and 
width («20 percent). 


CALCULATIONS 

Description of Nozzle Kinetic Program 

An IBM 7094 computer program for the solution of the one-dimensional flow of react- 
ing gases in a nozzle is described in references 15 and 16. Input information is in the 
form of tables of specific heat data, equilibrium constant data, nozzle area coordinates, 
and reaction rate data for a set of up to 15 reactions. The minimum number of reac- 
tions is equal to the total number of species minus the number of atom species. Starting 
conditions are initial values of pressure, temperature, velocity, and species concentra- 
tions. 

The program can be started at equilibrium or nonequilibrium conditions, and it will 
accept starting conditions either in the subsonic or supersonic portion of the nozzle. The 
problem of the singularity at the sonic point is handled by transforming the equation sys- 
tem in the transonic region to allow Mach number to be the independent variable rather 
than area. This results in slight alterations of the throat contour. 

Starting conditions for the kinetic program were computed from an equilibrium pro- 
gram (ref. 14). The equilibrium program computed combustor and nozzle conditions by 
starting with measured values of combustion air and fuel temperatures, pressure, and 
fuel-air ratio. The calculated equilibrium combustion temperature was not corrected 
for cooling losses to the combustion chamber and nozzle. Previous combustion temper- 
ature measurements (ref. 9) indicated that measured core temperatures were not af- 
fected by the cooled walls. Generally, starting conditions were taken from the equilib- 
rium calculations at subsonic area ratios of 1. 54 or 1. 04. In some cases where the ex- 
pansion remained closer to equilibrium, the kinetic program was started downstream of 
the throat at an area ratio of 1.01 to save computing time. 


Assumed Reaction Mechanism and Rate Constants 

Table I lists the standard set of reactions and rate constants used for comparison 
with the experimental results. The rates used for the hydrogen atom recombination re- 
actions (IV) and (V) were taken from reference 5 since they had shown good agreement with 
experimental temperature measurements when compared with the Bray freezing point 
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TABLE I. - RATES OF REACTIONS USED WITH KINETIC PROGRAM 


[units are cu cm, g mole, sec, and °K.] 


Reaction 

Standard rate 

Reference 

(I) H 2 + OH 32 H + H 2 0 

\ =(4.2x 10 12 )T 1 /2 e -10000/RT 

27 

(n) H + 0 2 * O + OH 

k n = (5.64xl0 13 )e' 1510 °/ RT 

27 

(III) O + H 2 - H + OH 

k m = ( 1 . 2xl0 13 )e' 9200/RT 

27 

(IV) H + H + M 52 H 2 + M 

hkpy = (3. 0 x10 18 )T -1 

5 

(V) OH + H + M 3t H 2 0 + M 

‘‘ky = (7. 5xl0 19 )T" 1 

5 

(VI) O + O + M at 0 2 + M 

k yI = (1.85xl0 17 )T -1 / 2 

27 

(VH) NO + M^N + O + M 

k yn = (1.06xl0 25 )T- 2 - 5 e- 154000/RT 

27 

(VIII) O + N 2 * N + NO 

k ym = (1.01xl0 12 )T 1//2 e" 7550 °/ RT 

27 

(DC) N + 0 2 « O + NO 

k K =(1.02xl0 11 )T 1 / 2 e- 620 °/ RT 

27 


a Rate variations were times 1/3 and 3. 

i. 

Rat e variations were times 1/10 and 100. 


analysis. Variations in the rates of reactions (1), (IV), and (V) were assumed to have the 
greatest effect on the overall kinetics (refs. 17 and 18). 


RESULTS AND DISCUSSION 
Hydroxyl Concentrations 

The experimental OH mole fractions, calculated by using equation (1), are given as 
a function of fuel-air equivalence ratio at nozzle supersonic area ratios of 1. 6, 2. 02, 

2. 49, and 3. 55 in figures 2(a) and (b). The stagnation pressures are 3. 6 and 1. 8 atmo- 
spheres, respectively. A dashed line was faired through the data for the Q^15 line. 

Discrepancies in the experimental data obtained at the higher area ratios were ob- 
served when different rotational lines were used in the transmission measurements. The 
higher energy lines (e. g. , Q^15) yielded the higher values of mole fraction. The differ- 
ences might be attributed to a combination of radial temperature profile and self- 
absorption, although an attempt at resolving the differences by assuming a temperature 
profile was not successful. 

The lower lines on figures 2(a) and (b) were calculated for equilibrium flow by using 
the IBM 7094 program of reference 14. Comparisons of the data with these equilibrium 
calculations indicate at least partial freezing of the OH mole fraction. It can be noted 
that the data show no significant effect of pressure and a much smaller dependence on 
equivalence ratio than the equilibrium curves. 
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OH mole fraction 



.4 .6 .8 1.0 1.2 .4 .6 .8 1.0 1.2 

Fuel-air equivalence ratio 


(a) Initial pressure, 3. 6 atmospheres; air temperature, 1710° to 1915° K. 


2x10' ^ 



Fuel-air equivalence ratio 

(b) Initial pressure, 1. 7 to 1. 9 atmospheres; air temperature, 1725° to 1865° K. 

Figure Z - Hydroxyl concentration by line absorption of ^TT transition; 0-0 vibrational band. Hydrogen-air combustion products. 


3000 


2800 1 


O Experimental 

□ Bray calculation with 



Figure 3. - Comparison of measured and calculated static tem- 
peratures. Variation of rate constant for H + OH + HoO + M; 
chamber pressure, 3. 6 atmospheres; equivalence ratio, 1.0. 


4xl0' 2 



Area ratio 


Figure 4. - Comparison of measured and calculated OH mole fraction. Variation 
of rate constant for H + OH + M * H 2 0+ M; chamber pressure, 3.6 atmo- 
spheres; equivalence ratio, 1.0. 



Comparisons with Kinetic Nozzle Calculations 


The experimental data are compared at a combustor pressure of 3. 6 atmospheres 
and at the stoichiometric fuel-air ratio with the kinetic calculations by using the standard 
rates and variations of the rate of reaction (V) (figs. 3 and 4). 

Figure 3 compares measured reversal temperatures along the axis of the nozzle 
(ref. 5) with temperatures calculated with the kinetic program by using variations of one- 
third, one, and three times the standard rate in reaction (V). The limiting cases of equi- 
librium and frozen flow and the calculations performed with the Bray freezing analysis 
are also indicated. The measured temperatures are seen to be about 70° K below the 
calculations that used the standard rates. A better comparison is indicated by ky times 
one-third. The Bray calculations that use the standard recombination rates compare 
quite well with the data, but they yield conservative values when compared to the com- 
plete kinetic calculation. 

Figure 4 compares the measured OH mole fraction with kinetic calculations. The 
data bars indicate the uncertainty in the measurement due to source intensity changes. 

If the differences introduced by using different rotational lines at the higher area ratios 
and the probable error in the f value and line shape corrections are considered, the 
measured OH mole fractions are consistent with the calculations that use values of ky 
between one-third and three times the standard rate. 

Figures 5 and 6 show comparisons of the data at 3. 6 atmospheres and at an equiva- 
lence ratio of 0. 6 with the kinetic calculations using the standard rates. The measured 
temperatures and OH mole fractions are in slightly better agreement than the data taken 
at <£> = 1.0. 

Comparisons of the data with calculations at a combustor pressure of 1. 8 atmo- 
spheres and at the stoichiometric fuel-air ratio are shown in figures 7 and 8. The mea- 
sured temperatures and Bray calculations are about 50° K below the kinetic calculations. 
Measured OH mole fractions are in only fair agreement with the calculations. 


Sensitivity of Calculations to Changes in Rate Constants 


The proximity of a reaction to equilibrium can be inferred by the parameter A when 


x j " k eq 


iIT c ‘ 
J 1 




(3) 


where is the concentration of specie i and v' and v" are the stoichiometric 
coefficients of reactants and products, respectively (in reaction j). The parameter A^ 
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OH mole fraction 


26001 — 



Figure 5. - Comparison of measured and calculated static tem- 
peratures. Standard reaction rates; chamber pressure, 

3. 6 atmospheres; equivalence ratio, 0. 6. 


2xl0' 2 



Area ratio 

Figure 6. - Comparison of measured and calculated OH mole fraction. 
Standard reaction rates; chamber pressure, 3. 6 atmospheres; equivalence 
ratio, 0.6. 


n 




Area ratio 

Figure 7. - Comparison of measured and calculated static tem- 
peratures. Standard reaction rates; chamber pressure, 

1. 8 atmospheres; equivalence ratio, 1.0. 


4xl0‘ 2 



Area ratio 


Figure 8. - Comparison of measured and calculated OH mole fraction Chamber 
pressure, 1. 8 atmospheres; equivalence ratio, 1.0. 
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TABLE H. - DEPARTURE FROM EQUILIBRIUM FOR 
VARIOUS REACTIONS WHEN STANDARD 
REACTION RATES ARE USED 


Area 

Lambda parameter 3. 

for first five reactions 

ratio 

x i 

A n 

X HI 

x rv 

V 

1.00 

1.00029 

1.00837 

1.02008 

1.27743 

1. 39532 

1.46 

1.00142 

1.01209 

1.03396 

14.7317 

17.8923 

1. 95 

1.00268 

.99934 

1.01522 

102. 259 

130. 555 

2.66 

1.00590 

.98234 

. 99445 

1 015. 50 

1 388. 24 

3. 19 

1.00935 

.97165 

.98513 

4 382. 15 

6 269. 06 

4. 16 

1.01917 

.95512 

. 98053 

42 336.7 

65 282.0 


0 th TT 

“For the j reaction, A. = K £ q | | Cj 

3 i 


is evaluated by the kinetic program and is listed in table n for the first five reactions as 
a function of area ratio. Because of the definition of the equilibrium constant, X = 1 at 
equilibrium. It is noted that reactions (I) to (III) remain near equilibrium throughout the 
area ratio range shown, whereas reactions (IV) and (V) depart from equilibrium upstream 
of the throat. 

Kinetic calculations were made to determine the effect of changing the rate constants 
of reactions (I) and (IV) (see table I, p. 7) . The calculations for temperature and OH mole 
fraction are shown in figures 9 and 10. No significant difference in temperature or OH 
mole fraction is introduced by changing the rate of the bimolecular reaction (I) by factors 
of one-third or three times the standard rate. This range is somewhat wider than the 
range of rate constants for reaction (I) found in the literature. The insensitivity to bi- 
molecular rates was expected in view of the assumption (ref. 18) that the bimolecular re- 
actions are in equilibrium with each other for the hydrogen- oxygen system. 

A decrease in the rate of the hydrogen atom recombination reaction by a factor of 10 
showed no significant change in the variables calculated by using the standard rates. An 
increase in this rate by a factor of 100, however, clearly allows this reaction to be con- 
trolling, and the calculated temperatures and mole fractions are shifted closer to equi- 
librium. Since the standard rate constant for reaction (V) was a factor of 25 greater 
than that of reaction (IV) , it is anticipated that significant effects would be noted if the 
rate of reaction (IV) were about equal or greater than reaction (V). 


Sensitivity of Calculations to Changes in Reaction Mechanism 


Fenimore and Jones (ref. 19) have inferred from some measurements in lean 


13 





Figure 9. - Effects of variation of rate constant data on static tempera- 
ture. Chamber pressure, 3. 6 atmospheres; equivalence ratio, 1.0. 


hydrogen- oxygen- argon flames that a recombination reaction involving HOg may be im- 
portant. The reaction scheme proposed is the three-body recombination reaction 


H + 0 2 + M * H0 2 + M 


AH = -47 kcal 


followed by the removal steps 


H0 2 + H - 2 OH 


AH = -38 kcal 


ho 2 + oh - o 2 + h 2 o 


AH = -38 kcal 


(X) 


(XI) 

(XII) 


14 



3xl0' 2 



Figure 10. - Effects of variation of rate constant data on OH mole fraction. 
Chamber pressure, 3. 6 atmospheres; equivalence ratio, 1.0. 


A rate for reaction (X) was inferred from low temperature data (refs. 19 to 21), an as- 
sumed T temperature dependence, and weighted third-body efficiencies (rate calcu- 
lated by Frank Belles, Lewis Research Center). The rate constant is given as 

k x = 7. 25 x 10 19 T -1 (cm 6 )(mole~ 2 )(sec" 1 ) 

The rate constants for reactions (XI) and (XII) were assumed to be 
k X i xn = 10 14 e~ 9090//RT (cm 3 ) (mole" ^ (sec" 1 ) 

Since the rate given in reference 19 is comparable with reaction (V), it would be 
likely to be of some importance at least in lean flames. 

Kinetic calculations were performed at equivalence ratios of 0. 6 and 1. 0 with the ad- 
dition of reactions (X) to (XII) to the reaction scheme. The rates used were taken from 
table I (p. 7) with ky times one and times one-third. The reduction in ky was used to 
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Figure 11. - Comparison of measured and calculated tempera- 
tures for a change in reaction mechanism; variation of rate 
constants for reactions V (H +OH +M^HoO +M) and X to 
XII(H+0 ? + /VUH0 2 +M, HOo + H-^2 OH, andH0 2 +0H* 
0 2 + H 2 0); chamber pressure, 3. 6 atmospheres; equivalence 
ratio, 0.6. 



Figure 12. - Comparison of measured and calculated OH mole fraction for a change 
in reaction mechanism; variation of rate constant for reactions V (H +OH +M 
H 2 0+M)andXtoXII<H+02+/VUH0 2 +M I H0 2 + H^20H, andH0 2 +0H^ 
0 2 + H 2 O); chamber pressure, 3. 6 atmospheres; equivalence ratio, 0.6. 


counter the effect on temperature of the additional recombination path since values of 
ky obtained from the literature would have included effects of k x if it has been impor- 
tant in the reaction scheme. 

The results of the calculations that used the alternate reaction mechanism are 
shown compared with the standard rates and the experimental data in figures 11 and 12 
for an equivalence ratio equal to 0. 6 and p c = 3.6 atmospheres. Results at an equiva- 
lence ratio of 1 . 0 and using ky times one showed negligible differences from the stan- 
dard set of reactions because of the low oxygen concentration. 
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At an equivalence ratio of 0. 6, the effect of including reactions (X) to (XU) more than 
compensated for the decrease in V Since the rates of reactions (X) to (XII) are uncer- 
tain, a further comparison was made with the standard set of rate constants and k x 
times one-tenth. The temperatures were not significantly different with or without reac- 
tion (X), but the OH mole fraction (fig. 12) showed a slight decrease. It appears that 
calculations which use some other combination of ky and k x might compare more 
favorably with the data at an equivalence ratio of 0. 6, but the uncertainties in the experi- 
mental data would not justify any conclusions as to their relative magnitudes. Further 
experimental determinations of the rates of reactions (V) and (X) to (XII) are needed for 
clarification of the recombination kinetics for lean mixtures. 


SUMMARY OF RESULTS 

Hydroxyl (OH) mole fractions were determined from spectral line absorption mea- 
surements in a supersonic nozzle expanding the combustion products of hydrogen and pre- 
heated air. Measurements were made as a function of fuel -air equivalence ratio (0. 6 to 

1. 2) at combustor pressures of 3. 6 and 1. 8 atmospheres at low supersonic area ratios 
(1. 6 to 3. 55). 

The measured OH mole fractions and previously reported static temperature mea- 
surements were compared to values of these parameters calculated by a one- dimensional 
computer program which included finite rate chemistry. The conclusions drawn from 
these comparisons are as follows: 

1. The OH radical concentration and temperature measurements compare reason- 
ably well with the kinetic computations that use a standard set of rate constants. The 
results of the OH concentration measurements are compatible with a rate constant for 
the reaction 


H + OH + M = H 2 0 + M 

within a factor of three higher or lower than the standard rate. The temperature mea- 
surements indicate a rate for the preceding reaction somewhat less than the standard rate. 

2. The insensitivity of the temperature and concentration to the bimolecular rates 
indicates that the partial equilibrium hypothesis is acceptable for the range of conditions 
investigated. 

3. Changes in OH mole fraction and temperature are sensitive to the recombination 
reactions 


H + H + MsH 2 + M 
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H + OH + M » H 2 0 + M 

The hydrogen recombination reaction only becomes important when its rate is com- 
parable to or larger than the rate of the water recombination reaction. 

4. For lean mixtures, the recombination reaction 

H + 0 2 + M * H0 2 + M 


followed by the removal steps 


H0 2 + H - 20H 
H0 2 + OH — 0 2 + H 2 0 

may be of importance since it offers an alternative recombination path. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 2, 1965. 
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APPENDIX A 


A/A* 



a 


b 

C 


c 

F 


/ 


h 

I 

J 

K 

K EQ 

k i 

k 


i 

m 

N 

P 

P* 


P 

Q R,V 


SYMBOLS 


nozzle supersonic area ratio 

relative line strength 

absorptance 

broadening parameter 

half width at half intensity, cm 

concentration, moles/cu cm 

speed of light, cm/sec 

^(2J ,f + 1)/Ag, function 

ratio of number of dispersion 
electrons to number of ab- 
sorbers 

Planck constant 
relative intensity 
rotational quantum number 



equilibrium constant 

rate constant, j = I, II • • • XII 

Boltzmann constant 

path length, cm 

atomic mass 

_3 

number density, cm 

absorption coefficient 

absorption coefficient for a 
Doppler shaped line 

pressure, atm 

rotational-vibrational partition 
function 


R gas constant 

T temperature, °K 

Tj«j»» correction factor to line 
strengths for vibration- 
rotation interaction 

y variable of integration 

a ratio of source line to absorber 

line widths 

X defined by eq. (3) 

v stoichiometric coefficient 

(p fuel-air equivalence ratio 

w wave number, cm - * 

Subscripts: 

A absorption line 

C collision 

c chamber 

D Doppler 

E emission line 

i species 

j reaction 

K rotational line designation 

N natural 

a) wave number 

ci) Q wave number at line center 

0 initial 

Superscripts: 

’ reactant 

,f products 


19 



APPENDIX B 


EVALUATION OF LINE SHAPE CORRECTIONS 

Equation (1) (p. 4) is exact when P T represents the absorption coefficient at the 
center of a Doppler shaped line. Since the OH absorption lines are partially collision 
broadened and the source emission lines have finite width, corrections must be made to 
the absorption coefficient in equation (1) so that the integral over the absorption line rep- 
resents the value for the actual line shape. 

The absorptance is related to the absorption coefficient at the line center by (ref. 22, 
pp. 122 and 169) 



and 


_ Aa) E _ b D(E) 
a “A ( b C +b D) A 


(B3) 


with bp and b^-, as the Doppler and collision half widths at half intensity. 

The collision half width of the absorption line can be written in terms of the broad- 
ening parameter a where 


a = (In 2)1/2 (B4) 

b D 

When the natural half width b^ is neglected in comparison with b^,, equation (B3) be- 
comes 
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a . - b D(E) 


1 


3 D(A) 1 + a(ln 2) 


• 1/2 


(B5) 


For optically thin emission and absorption lines, the ratio of Doppler half widths is 
equal to the square root of the absolute temperature ratio and 


a = 



T A 1 + a(ln 2) -1 / 2 


(B6) 


Equation (Bl) has been evaluated for Doppler shaped absorption lines in appendix IV 
of reference 22, where the absorption term is 1 - expj^-P'f exp^-co 2 ^ . 

It can also be evaluated from the series expansion 

,2 


A.. = 


P'f 


(p'lV 


/ 2\^ 2 ( 2\^ 2 

(l + a A ) 2! (l + 2cT) 


+ (-l) 


n 


(P'l) n 


n! ^1 + na 


2)V2 


(B7) 


For combined Doppler and collision broadening, the absorption coefficient is related to 
the coefficient at the line center by (ref. 22) 


-y = ^ f [a 2 + (w - y) 2 ] [exp(-y 2 )] dy 

t/-oo 


(B8) 


where 


P dco 


, 1/2 


P’ = 


•/-oo 

/ me 

CO 

\2tt/T 


(B9) 


is the absorption coefficient at the center of a Doppler broadened line. 

For small values of a equation (B8) can be written as an expansion around the 
Doppler contour (ref. 23): 


= exp^-co 2 ^ - 2a" 1//2 |/. - 2w 1 F(u> 1 )J + a 2 (l - 2o> 2 )exp (-u> 2 ) 

-2a 3 - (l - o> 2 )- 20)^1 - 1 cu^F^j) 


(BIO) 
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In order to determine the proper value of P* to use in equation (1), an approximate 
approach was used here in which the corrections for finite emitter line width and ab- 
sorber line shape were evaluated separately. As a first approximation, the absorption 
line was assumed to be Doppler broadened only, and the right side of equation (Bl) eval- 
uated from the approximate relation (B7) and the tabular values of appendix IV in refer- 
ence 22. Values of £ were plotted as functions of 


I 1 - A 


a 


for various values of a . The value of the absorption coefficient is obtained from the 

plot by using experimental values of Iq/I, a, and the path length Si. Equation (BIO) is 

evaluated at the line center where u>, = 0 and where P v\ the ratio of the absorp- 

o/ 

tion coefficient for the actual line shape to the coefficient for a Doppler shaped line, is 

obtained. When the values determined for P /t and P /p' are used, the appropriate 

o %/ 

value for P* can then be substituted in equation (1). The total correction increased the 
value of the absorption coefficient by about 20 percent at a chamber pressure of 3. 6 atmo 
spheres and 17 percent at a pressure of 1. 8 atmospheres. 

The error in using the Doppler form of equation (Bl) was evaluated for the experi- 
mental data having the largest values of the broadening parameter a by numerical inte- 
gration using equations (Bl) and (BIO). The approximate correction was found to yield 
values of the absorption coefficient which were a maximum of 10 percent low, and most 
of the data had errors in the absorption coefficient which were well under this value. In 
view of much larger uncertainties in the parameters a and a, use of the more exact ex 
pression was not justified. 

The broadening parameter a was evaluated from (ref. 24) 


a = 335 E- (Bll) 

T 

and from measured values of the static pressure and temperature. Equation (Bll) is an 
empirical relation derived from the pressure and temperature dependence of the colli- 
sion and Doppler widths. 

Carrington (ref. 25) has evaluated the parameter a by using the curve of growth 
and has reported a = 0. 05±0. 05. At the temperature of his flame, 1 equation (Bll) gives 
a value of a = 0. 117. Kaskan (ref. 26) has given a = 450 £_ as best fitting his experi- 

^Carringtons* measured rotational temperature of 2660° K has been corrected to 
2866° K by Learner (ref. 11) for vibration- rotation interaction. 
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mental measurements at temperatures around 1500° K. Considering the uncertainty in 
the values of references 25 and 26, equation (Bll) is consistent with either. 

The value of a is calculated from equation (B6) by using the measured reversal 
temperatures for and a source translational temperature taken somewhat arbitrarily 
as T E = 500° K. The source rotational distribution could not be characterized by a sin- 
gle temperature, but the distribution curve appeared to have an asymptotic slope for ro- 
tational quantum numbers less than or equal to K = 3, which gave the above value. A 
value of T E = 450° K has been reported for a similar discharge source (ref. 13). 
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